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Abstract 

We present two recent applications of lattice-gas modeling techniques to elec- 
trochemical adsorption on catalytically active metal substrates: urea on Pt(100) 
and (bi)sulfate on Rh(lll). Both systems involve the specific adsorption of small 
molecules or ions on well-characterized single- crystal electrodes, and they provide a 
particularly good fit between the adsorbate geometry and the substrate structure. 
The close geometric fit facilitates the formation of ordered submonolayer adsor- 
bate phases in a range of electrode potential positive of the range in which an 
adsorbed monolayer of hydrogen is stable. In both systems the ordered-phase re- 
gion is separated from the adsorbed-hydrogen region by a phase transition, signified 
in cyclic voltammograms by a sharp current peak. Based on data from in situ ra- 
diochemical surface concentration measurements, cyclic voltammetry, and scanning 
tunneling microscopy, and ex situ Auger electron spectroscopy and low-energy elec- 
tron diffraction, we have developed specific lattice-gas models for the two systems. 
These models were studied by group-theoretical ground-state calculations and nu- 
merical Monte Carlo simulations, and effective lattice-gas interaction parameters 
were determined so as to provide agreement with the experimental results. 
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1 Introduction 



Both single-component and multi-component adsorption processes are strongly 
influenced by lateral adsorbate-adsorbate interactions. A theoretical approach 
to the investigation of these effects is the formulation of lattice-gas models and 
their solution by analytical and numerical techniques. This method has previ- 
ously been successfully applied to phase transitions at solid- vacuum and solid- 
gas interfaces [1]. Provided care is taken to consider its limitations, the lattice- 
gas modeling approach is also useful to understand the statistical-mechanical 
aspects of multicomponent adsorption at solid-liquid interfaces, in particular 
the complicated adsorbate phase diagrams observed near the potential of zero 
charge [2-13]. 

In this paper we discuss the application of computational lattice-gas modeling 
to two different experimental systems involving adsorption on well-characterized 
single-crystal electrodes: urea on Pt(100) [14-19] and (bi)sulfate on Rh(lll) 
[20-22]. Both systems exhibit a dramatic peak sharpening in the cyclic voltam- 
mogram (CV) from several hundred mV to on the order of 10 mV when the 
adsorbate species (urea or (bi)sulfate) is added to the supporting perchloric- 
acid electrolyte [14,20]. We associate this effect with a phase transition in the 
layer of contact adsorbed particles. This transition involves the replacement 
of a monolayer of adsorbed hydrogen on the negative-potential side of the CV 
peak by an ordered submonolayer of urea or (bi)sulfate, respectively, on the 
positive-potential side. 

The peak-sharpening phenomenon is much weaker or absent when the same 
substances are electrosorbed onto other crystal planes of the same metals 
[14,20]. This high specificity with respect to the adsorbent surface structure 
indicates that the effect depends crucially on the geometric fit between (at least 
one of) the adsorbate species and the surface. Guided by this observation, we 
have developed specific lattice-gas models with three adsorbed species: wa- 
ter, hydrogen, and either urea or (bi)sulfate, depending on the system. These 
models were studied by group-theoretical ground-state calculations [23-25] 
and numerical Monte Carlo (MC) simulations [26]. Based on the results of 
these calculations, we adjusted the effective lateral adsorbate-adsorbate inter- 
actions in the models to fit experimentally measured quantities. In addition 
to the systems discussed here, similar modeling strategies, differing mainly 
in the means by which thermodynamic information was extracted from the 
lattice-gas Hamiltonian, have been applied to other systems involving local- 
ized adsorption of two species of small molecules or ions [3-12]. For the sys- 
tems discussed here, the experiments with which agreement was sought include 
cyclic voltammetry, in situ radiochemical surface concentration measurements 
(RCM), and ex situ Auger electron spectroscopy (AES) and low-energy elec- 
tron diffraction (LEED). 
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Our lattice-gas models are defined through a generalization of the standard 
three-state lattice-gas Hamiltonian (energy function) used, e.g., in Refs. [3-6]: 
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Here cf £{0,1 } is the local occupation variable for species X [X=A (adsorbate) 
or H (hydrogen)], and the third adsorption state ("empty" or "solvated") cor- 
responds to cf =cf=0. The sums Y^ifj) an d J2i run over all nth-neighbor bonds 

and over all adsorption sites, respectively, $xy denotes the effective XY pair 
interaction through an nth-neighbor bond, and J2 n runs over the interaction 
ranges. The term 7i 3 contains three-particle [27] and possibly multi-particle 
interactions. Both the interaction ranges and the absence or presence of multi- 
particle interactions depend on the specific system. The change in electro- 
chemical potential when one X particle is removed from the bulk solution and 
adsorbed on the surface is —fix- The total number of surface unit cells is iV, 
and ©x = iV -1 J2i 4 i s the surface coverage by species X. The sign conven- 
tion is such that $xy>0 denotes an effective attraction, and j2x>0 denotes 
a tendency for adsorption in the absence of lateral interactions. Straightfor- 
ward modifications of the above form are necessary if the adsorption sites for 
the two species are different, as they are in the two cases studied here. We 
emphasize that the interactions in Eq. (f ) are effective interactions mediated 
through a variety of channels, including electrons, surface phonons, and the 
fluid [5], and that they therefore in principle may depend on temperature and 
electrode potential. 

The electrochemical potentials in Eq. (f) are (in the weak-solution approxi- 
mation) related to the bulk concentrations [X] and the electrode potential E 
as 

Ji x = fil + RThi^-z x FE, (2) 



where R is the molar gas constant, T is the absolute temperature, F is Fara- 
day's constant, and the effective electrovalence of X is zx- The quantities 
superscripted with a are reference values which contain the local binding 
energies to the surface. 

For a given set of effective interactions, the coverages ©x as functions of the 
electrochemical potentials and temperature are calculated numerically by MC 
simulation. In the absence of diffusion and double-layer effects and in the limit 
that the potential sweep rate dE/dt— s-0, the CV current per unit cell of the 
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surface (oxidation currents considered positive) is related to the lattice-gas 
response functions d®x/djj>Y as [16-18] 



2 <90 A . <90 H 2 <90 H | dE 



where e is the elementary charge unit. 

Once a lattice-gas Hamiltonian has been specified, the corresponding equilib- 
rium thermodynamic properties can be studied by a number of analytical and 
numerical methods, depending on the quantities of interest and the complexity 
of the Hamiltonian. In addition to group-theoretical ground-state calculations 
to determine the domains of stability of the various possible phases [23-25], 
finite-temperature properties have been obtained by mean-field approxima- 
tions [7] (although these can be unreliable for low-dimensional systems with 
short-range interactions [28,29]), Pade-approximant methods based on liquid 
theory [9-12], numerical transfer-matrix calculations [3-6], and Monte Carlo 
simulations [4,16-18]. 



2 Urea on Pt(100) 



In addition to the aforementioned, surface-specific narrowing of the CV peak 
upon the addition of urea to the supporting electrolyte, the experimental 
observations to which we have attempted to fit the model are as follows. (For 
experimental details, see Ref. [18].) 

1. The urea coverage ©u, measured in situ by RCM, changes over a potential 
range of approximately 20 mV around the CV peak position from near zero on 
the negative side to approximately 1/4 molecules per Pt(100) unit cell (1/4 
monolayers, or ML) on the positive side (0.26±0.04ML at — 60 mV versus 
Ag/AgCl), as shown in Fig. lb. Corresponding CV curves are shown in Fig. la. 

2. Ex situ AES studies of the electrode emersed at — 60 mV are consistent with 
the RCM result (0 U= O.24±O.O3 ML). 

3. Ex situ LEED studies of the electrode emersed at — 60 mV show an ordered 
c(2x4) adsorbate structure, consistent with an ideal coverage of 1/4 ML. Upon 
emersion on the negative side of the CV peak, only an unreconstructed (lxl) 
surface (presumably prevented from reconstructing by trace impurities) was 
found. 

4. The CV peak potentials were measured for [U] between 0.3 and 3.0 mM. In 
Fig. 2a are shown the corresponding electrochemical potentials ( X and dotted 
line), obtained from Eq. (2) with the same parameters used to fit Fig. 1. 

5. In addition to the above experiments, which were all performed at room 
temperature, cyclic voltammetry was performed at several fixed temperatures 
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Fig. 1. Urea on Pt(100) in 0.1 M HCIO4 at room temperature, (a) Experimen- 
tal (dashed curves) and simulated (O and solid curve) normalized CV currents, 
i/(dE/dt) in elementary charges per mV per Pt(100) unit cell, at 1.0 mM bulk 
urea. The two dashed curves are representative negative-going voltammograms, and 
their differences indicate the experimental uncertainty, (b) Urea coverage measured 
by RCM at 0.5 mM bulk urea (x with error bars) and simulated by MC (O and 
solid curve), shown versus electrode potential. Also shown as a dotted curve is the 
simulation result for 1.0 mM bulk urea, which corresponds directly to the simu- 
lated CV profile shown in (a). The effective interactions used in this work are (in 
kJ/mol) $$ = -2.0, $$ = -8.0, $$=-4,0, $$ = -13.0, $$ = -10.0, $$=-5.9, 

$$, = -0.5, $$ = -2.5, $$ = -3.0, $$ =+0.25, $$ = -2.0, the effective elec- 
trovalences were taken as zn = + l and z\j = — l, and the room-temperature simula- 
tions were performed at 28°C. The results shown are for a system of 32x32 Pt(100) 
unit cells. The statistical uncertainty in the MC results is everywhere much smaller 
than the symbol size. 

between 0°C and 40° C. 

The lattice-gas model developed to account for these observations is based 
on the assumption that urea [CO(NH 2 )2] coordinates the platinum surface 
through its nitrogen atoms (or NH 2 groups), with the C=0 group pointing 
away from the surface. Since the unstrained N-N distance in urea matches 
the lattice constant of the square Pt(100) surface quite well (2.33 A [30] ver- 
sus 2.77 A [31]), we assume that urea occupies two adsorption sites on the 
square Pt(100) lattice. Integration of the CV curves indicates that the hydro- 
gen saturation coverage in the negative-potential region corresponds to one 
elementary charge per Pt(100) unit cell, and that most of the surface hydro- 
gen is desorbed in the same potential range where urea becomes adsorbed. 
Therefore, we assumed [15] that hydrogen adsorbs in the same on-top posi- 
tions as the urea nitrogen atoms. This assumption was recently confirmed by 
visible-infrared sum generation spectroscopy [32]. The resulting model [15] is 
a dimer-monomer model in which hydrogen is adsorbed at the nodes and urea 
on the bonds of a square lattice representing the Pt(100) surface. Simulta- 
neous occupation by urea of two bonds that share a node is excluded, as is 
occupation by hydrogen of a node adjacent to a bond occupied by urea. In or- 
der to stabilize the observed c(2x4) phase, effective interactions are included 
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Fig. 2. Phase diagrams for urea on Pt(100) in 0.1 M HCIO4, shown in the (fiv^fiu) 
plane, (a) The zero-temperature phase boundary between the c(2x4) phase with 
©u = l/4, indicated as c(2x4)U in the figure, and the (lxl) phase with 0^ = 1, 
indicated as (lxl)H, is shown as a solid line together with the electrochemical po- 
tentials corresponding to experimental (x connected by dotted lines) and simulated 
( + connected by dashed lines) CV peak positions. All simulation parameters are 
as in Fig. 1. The solid arrows represent positive-going E scans from — 106 mV to 
— 56 mV. From left to right they represent [U]=2mM, ImM, and 0.5 mM. The ex- 
perimental peak positions are averages over positive-going and negative-going scans 
with |di?/d/| = 50mV/s, but the average positions depend little on |di?/d/|. (b) A 
full ground-state diagram for the interactions used, showing all the phases present. 
The phase indicated in the figure by c(2x4)UU has ©u = l/2, (lx3)U has u = l/3, 
(v^Xv^)!! has ©h = 1/2, and (1x1)0 is the empty lattice. The phase regions out- 
side the dotted box, which corresponds to panel (a), are likely not experimentally 
accessible. 



through eighth-nearest neighbors [16-19], as shown in Fig. 5 of Ref. [18]. The 
model is illustrated in Fig. 3. 



The effective lattice-gas interactions (given in the caption of Fig. 1) were de- 
termined by a ground-state calculation followed by numerical MC simulation. 
Since the number of model parameters is large, the numerical calculations 
are time consuming, and the experimental data concern a number of different 
quantities, parameter estimation by a formal optimization procedure was not 
a practical alternative. (This contrasts with the simpler situations discussed 
in Refs. [5,29,33], where a small number of lattice-gas parameters could be 
determined by a formal least-squares procedure to fit extensive experimen- 
tal results for a single thermodynamic quantity.) The model parameters were 
therefore varied "by hand", taking into consideration both the various exper- 
imental results and available chemical and physical background information, 
until acceptable agreement was obtained with room-temperature experimental 
results. In particular, we sought agreement between the simulated CV profiles 
and those obtained experimentally at the lowest practical potential scan rate 
of 0.5 mV/s, while at the same time ensuring that the simulated ordered phase 
on the positive side of the CV peak corresponded to the observed c(2x4) struc- 
ture with 1/4 MF urea coverage. Whereas this procedure may be criticized as 
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Fig. 3. Left panel: Illustration of the microscopic model, showing urea molecules 
with their NH2 groups (dark gray) on the square Pt(100) lattice sites and their C = 
groups (C black and lighter gray) pointing away from the surface. The hydrogen 
atoms are shown as single light gray spheres. Right panel: A typical equilibrium 
configuration generated by MC in the ordered c(2x4) phase region at —60 mV. Urea 
molecules are shown as black rectangles on the lattice bonds and hydrogen as black 
circles at the nodes. Locally disordered fluctuations in the ordered background are 
seen, as well as a low density of hydrogen on vacant sites between the urea molecules. 
This figure is adapted from the video, Ref. [19]. 

unsystematic, we believe it makes the maximum use of all the available exper- 
imental and theoretical information. A detailed discussion of the roles of the 
different effective interactions is given in Ref. [18]. 

The numerical simulations, which used systems with up to 32x32 square- 
lattice unit cells, were performed on a cluster of IBM RS6000 workstations. 
We used a heat-bath MC algorithm [26] with updates of clusters consisting 
of five nearest-neighbor nodes arranged in a cross, plus their four connecting 
bonds. After symmetry reductions these clusters have 64 different configura- 
tions, and the corresponding code is rather slow in terms of machine time per 
MC step. However, the additional transitions allowed by these clusters, relative 
to minimal clusters consisting of two sites and their connecting bond, include 
"diffusion-like" moves in which the urea molecules can go from one bond to 
another and the hydrogen atoms from one site to another, without changing 
the local coverages within the cluster. These moves significantly reduce the 
free-energy barriers that must be surmounted in order to locally minimize the 
adsorbate energy, and consequently dramatically reduce the number of MC 
steps per site (MCSS) necessary for the system to reach thermodynamic equi- 
librium. This speedup is particularly evident deep in the c(2x4) phase, where 
an algorithm based on updates of minimal clusters proved unable to equili- 
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(a) 



(b) 




Fig. 4. Simulated "LEED patterns" in the Brillouin zone surrounding a Bragg peak. 
The simulated scattering intensity is shown on a logarithmic scale versus the x- and 
components of the dimensionless scattering vector, q x and q y . Contour plots of 
the intensity are shown in the (q x ,q y ) plane. Both panels are for ImM urea, (a): 
E = — 0.72 mV. (b): E = — 0.81 mV. See discussion in the text. 



brate the system, even after several million MCSS. Thermodynamic quantities 
were calculated as averages over each MC run, with samples taken every fifty 
MCSS. 

Near the CV peak at —91 mV, the simulated height and width of which agree 
well with the room-temperature experimental results as shown in Fig. la, no 
finite-size effects were observed for system sizes above 8x8. (For a discussion 
of temperature effects, see below.) This indicates that this peak in the model 
is not directly due to a phase transition [34], but probably results from strong, 
local coverage fluctuations in a disordered phase near a tricritical point. At 
the slightly more positive potential of — 75 mV a second-order, order-disorder 
transition (probably in the universality class of the XY model with cubic 
anisotropy) occurs, which is indicated by a sharp peak in the response func- 
tion for the c(2x4) order parameter. That this peak indeed signifies a phase 
transition, is indicated by its finite-size scaling behavior [34]. Since it does not 
involve large changes in the coverages, in the room-temperature simulation 
this order-disorder transition only shows up in the CV as a small shoulder on 
the positive side of the peak. 

Simulated "LEED patterns" were obtained as the squared Fourier transform 
of the adsorbed urea configurations. These were obtained by the Fast Fourier 
Transform algorithm (our program was based on the subroutine fourn of 
Ref. [35]) and averaged in the same way as the thermodynamic quantities. 
In Fig. 4 we show results for ImM urea at — 72mV, which is in the c(2x4) 
region on the positive side of the order-disorder transition (Fig. 4a), and at 
—81 mV, which lies between the order-disorder transition and the CV peak 
(Fig. 4b). The triple intensity maxima at the edges of the Brillouin zone, which 
are characteristic of c(2x4) order, are clearly seen. In the ordered phase they 
are comparable in intensity with the peak at the zone center, the intensity of 
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Fig. 5. Temperature effects in experiments and model, (a) The measured temper- 
ature dependence of the full width at half maximum of the CV peak (FWHM) 
during positive-going (O) and negative-going (□) scans, and the peak separation 
(x). For clarity, the FWHM data points are offset slightly to the right (O) and left 
(□) from the actual temperature values. Along the upper edge we give the number 
of independent scans at each temperature, and the error bars correspond to the em- 
pirical standard deviation, (b) Simulated CV profiles for a 32x32 system at 28°C 
(O, identical to Fig. la) and 10°C (x and solid line). All model parameters are as 
in Fig. 1. 



which is the square of the urea coverage. In contrast, the c(2x4) maxima in 
Fig. 4b are much weaker than the peak at the zone center. This indicates that 
the long-range order is almost totally absent at this potential, even though 
the urea coverage is near that of the ordered phase. 

Exploratory cyclic voltammetry experiments have been carried out between 
0°C and 40° C. The temperature was adjusted by immersing the electrochemi- 
cal cell in a thermostat controlled water bath. A number of independent poten- 
tial cycles were carried out at each temperature. In Fig. 5a we show the results 
for the full width at half maximum (FWHM) and the potential separation be- 
tween the positive-going and the negative-going peak at |d.E/dt| = l.OmV/s. 
Only temperatures for which at least four potential cycles were performed are 
included. The data shown are averages over the number of cycles indicated 
along the upper edge of the figure, and the error bars represent one empirical 
standard deviation. No statistically significant temperature dependence was 
found for the FWHM. The same was the case for the maximum currents and 
the average peak positions, neither of which are shown here. By contrast, the 
peak separation at fixed |d.E/dt| = l.OmV/s was reduced approximately by 
a factor of six from the lowest to the highest temperature shown. Whereas 
this reduction is highly statistically significant, the temperature range is too 
small to distinguish reliably between, e.g., a linear decrease and an Arrhenius- 
type behavior on the basis of these experiments alone. However, the reduction 
with increasing T of the peak separation at constant scan rate indicates the 
presence of significant temperature-dependent hysteresis, even at this low scan 
rate. 
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For comparison with these experimental results, we present in Fig. 5b sim- 
ulated CV profiles at 10°C and 28°C (the latter identical to the one shown 
in Fig. la). Notice that the shape and height of the non-scaling peak which 
dominates the simulated CV current at 28° C is hardly changed at all when the 
temperature is lowered to 10°C. In contrast, the scaling peak corresponding 
to the order-disorder phase transition, which at 28° C appears only as a small 
shoulder near the foot of the main peak, at 10° C appears as a sharp peak 
with an integral corresponding to approximately 15% of the total transported 
charge. 

The simulated CV profile at this temperature does not agree well with the 
experimentally observed profiles, which do not differ materially from those of 
the room-temperature voltammograms shown in Fig. la. The main difference 
between the simulated profiles at the two temperatures lies in the sharp spike 
associated with the order-disorder phase transition. We believe this discrep- 
ancy between the simulated and experimental CV profiles at low temperatures 
can be explained by the temperature-dependent kinetic effects evidenced by 
the experimental peak separations shown in Fig. 5a. Such effects could in prin- 
ciple be due to bulk diffusion, as well as to surface kinetics including diffusion 
and ordering. For the bulk concentration range used in this study, bulk dif- 
fusion should be unimportant, and we therefore believe the hysteresis is due 
to processes on the surface. It is difficult to quantify the influence of these 
kinetic effects on the CV peak shapes and positions without either employing 
highly phenomenological approximations or carrying out extensive nonequi- 
librium simulations. Qualitatively, however, one can point out that whereas 
the characteristic times associated with the wide peak, which reflects only lo- 
cal coverage fluctuations on the surface, should be relatively short (estimated 
from the experimental peak separations as 1-10 s), the timescale associated 
with the narrow ordering peak should be considerably longer. Two distinct 
physical phenomena contribute to this effect. The first is the critical slowing- 
down associated with the continuous order-disorder phase transition [26]. The 
second is a "gridlock" phenomenon [36] which prevents a disordered urea layer 
with a coverage near 1/4 ML from easily adjusting to the ordered equilibrium 
configuration, due to the large free-energy barriers that prevent desorption 
that could momentarily ease the crowding. This effect is closely related to the 
dramatic slowing-down suffered in the ordered c(2x4) phase region by the MC 
algorithms that do not allow lateral diffusive moves. As a result, the narrow 
peak in the simulated "equilibrium" CV profile should be much more seri- 
ously affected by the surface kinetics, to the extent of becoming smeared into 
an almost undetectable shoulder on the side of the broader peak in the exper- 
imental voltammogram. (Recent reviews of related issues of metastability and 
hysteresis in a variety of fields are found in Refs. [37,38].) 

In conclusion, concerning the discussion in the previous paragraph, we believe 
that the observed kinetic effects are too large to allow us to determine whether 
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or not the sharp ordering peak predicted by the lattice-gas model at the lower 
temperatures is present in the equilibrium behavior of the experimental sys- 
tem. We can suggest several possible methods by which one could resolve this 
question. First, voltammetry could be performed at lower scan rates in order 
to attempt to resolve any sharp ordering peak. However, the required scan 
rates may well be too low to be practical [39]. Second, one could relegate the 
CV data to a secondary position in the hierarchy of experimental results to 
be fitted, concentrating instead on RCM coverage measurements, which are 
obtained much closer to equilibrium conditions. This is the approach taken in 
our study of (bi)sulfate adsorption on Rh(lll), discussed in Sec. 3. However, 
the uncertainties in the RCM coverage measurements (5-10%) and electrode 
potential determination (ss ±2mV) makes detection by this method of a nar- 
row peak representing a small total charge very difficult. Third, one could 
attempt to account for the effects of the surface kinetics on the CV peak 
shapes and separations, either analytically, as it was recently done in a study 
of underpotential deposition of copper on Au(lll) in sulfuric acid [11], or 
by kinetic MC simulations. However, the MC algorithms used in the present 
study were designed to equilibrate the system in a reasonably short amount 
of computer time, rather than to faithfully reproduce the actual ordering ki- 
netics. We therefore have not attempted to systematically study the influence 
of the lattice-gas kinetics on the observed voltammograms at finite scan rates. 
This problem will be considered in the future. 



3 (Bi)sulfate on Rh(lll) 



In this section we report some preliminary results of our ongoing study of 
(bi)sulfate electrosorption on Rh(lll) single crystals. Details will be reported 
elsewhere [21,22]. In addition to the surface-specific narrowing of the CV peak 
upon the addition of sulfuric acid to the supporting electrolyte, the experi- 
mental observations to which we have attempted to fit a lattice-gas model are 
as follows. 

1. The (bi)sulfate coverage 0s, measured in situ by RCM, changes over a po- 
tential range of approximately 15 mV around the CV peak position from near 
zero on the negative side to approximately 1/5 molecules per Rh(lll) unit cell 
(1/5 ML) on the positive side (0.21±0.01 ML at -140 mV versus Ag/AgCl), 
as shown in Fig. 6. 

2. In situ scanning tunneling microscopy (STM) studies of (bi)sulfate on 
Au(lll) [40-42] and Rh(lll) [42] indicate an ordered (\/3x\/7) sulfate sub- 
monolayer at potentials sufficiently far to the positive side of the CV peak. 
The most likely (bi)sulfate coverage deduced from these studies is approxi- 
mately 1/5 ML, possibly with additional coadsorbed water or hydronium ions 
[41]. 
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Fig. 6. (Bi)sulfate coverages in 0.1 M HCIO4 supporting electrolyte, as measured 
by the in situ RCM technique for bulk sulfuric acid concentrations 0.1 mM (x), 
0.3 mM (O), 1.0 mM (□), and 3.0 mM ( + ). Together with the experimental data 
are shown preliminary MC simulation results for 0.3 mM (solid curve) and 1.0 mM 
(dashed curve), using the simplified lattice-gas model described in the text. 



3. Ex situ AES studies of the electrode emersed at — 140 mV are consistent 
with the RCM result, yielding values of 0s near 1/5 ML [43]. 

4. Ex situ LEED studies of the electrode emersed at — 140 mV show an or- 
dered (y/3xy/3) structure [43]. The ideal coverage of this structure (1/3 ML) 
is not consistent with the in situ RCM and STM and the ex situ AES results. 
However, it is possible that as the (bi)sulfate adsorbate system becomes an ef- 
fectively closed system upon emersion, the change of environment may change 
the effective interactions (possibly through desorption of strongly coadsorbed 
water or hydrogen) to cause condensation into (y/3xy/3) islands coexisting 
with patches of essentially empty surface, to give a total coverage near 1 /5 ML. 
This explanation has recently been suggested for analogous observations for 
(bi)sulfate on Au(lll) [44]. 

5. The CV peak potentials were measured for bulk sulfuric acid concentrations 
between 0.1 and 10 mM. 



The lattice-gas model developed to account for these observations is based on 
the assumption that the sulfate (SO^ - ) or bisulfate (HSO4 ) coordinates the 
Rh(lll) surface through three of its oxygen atoms, with the fourth S-0 bond 
pointing away from the surface, as is also the most likely adsorption geometry 
on Au(lll) [8-12,41]. This adsorption geometry gives the sulfate a "footprint" 
in the shape of an approximately equilateral triangle with an 0-0 distance of 
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2.4 A [45], reasonably matching the lattice constant for the triangular Rh(lll) 
unit cell, 2.69 A [31]. This high symmetry is somewhat distorted for bisulfate. 
Based on spectroscopic information, it is likely that the (bi)sulfate coordi- 
nates the rhodium atoms directly through three oxygen atoms, thus bringing 
the sulfur directly above the three-fold hollow sites [46]. Integration of the 
CV curves indicates that the hydrogen saturation coverage in the negative- 
potential region corresponds to one elementary charge per Rh(lll) unit cell, 
and that most of the surface hydrogen is desorbed in the same potential range 
where (bi)sulfate becomes adsorbed. As for platinum [32], we expect that the 
hydrogen is adsorbed at on-top sites. In the resulting model the (bi)sulfate is 
modeled as hard hexagons [8-12,47], and the hydrogen as point particles, and 
finite effective pair interactions out to fourth-nearest neighbors are included 
to stabilize the (y/3xy/7) phase. In addition, three-particle interactions [27] 
are included to discourage the formation of the (y/3xy/3) phase. 

The effective lattice-gas interactions can be estimated as described for urea on 
Pt(lll) in Sec. 2, and simulations were performed on the workstation cluster 
described in Sec. 2 for systems up to 30x60 triangular-lattice unit cells. For the 
same reasons of computational speed discussed in Sec. 2, we used a heat-bath 
algorithm with updates of clusters consisting of two nearest-neighbor sites, 
which have a total of six possible adsorption states. Preliminary simulation 
results for 0s in a simplified model, in which the hydrogen and the sulfate 
are competitively adsorbed on the same sites, are shown in Fig. 6, together 
with the experimental RCM results. These simulation results agree quite rea- 
sonably with the (bi)sulfate saturation coverage, as well as with the general 
shape of the experimental coverage curves in the transition potential region 
and with the dependence of the transition potential on the bulk concentra- 
tion of sulfuric acid. The "jaggedness" of the simulated coverage curves in the 
0.17-0.20 ML region is a consequence of the relatively short simulation runs 
used (2.4 xlO 5 MCSS per point). Further simulations are in progress. A rep- 
resentative equilibrium configuration for the simplified model in the ordered 
(y/3xy/7) phase region is shown in Fig. 7. 



4 Discussion 

In this paper we have discussed two particular applications of lattice-gas mod- 
els with effective lateral interactions to electrochemical adsorption of small 
molecules and ions on well-characterized single-crystal surfaces. Detailed ex- 
perimental information is available for both systems studied. A characteristic 
aspect of these systems is the good geometric fit between the adsorbent surface 
lattice and the main adsorbate. This allows the formation of ordered adsor- 
bate phases commensurate with the substrate, which can be observed both 
by in situ atomic-scale microscopies and by ex situ scattering techniques. The 
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Fig. 7. A typical equilibrium configuration for our simplified model for (bi)sulfate 
adsorption on Rh(lll), generated by MC in the ordered (\/3X\/7) phase region at 
— 150 mV. (Bi)sulfate ions are shown as O? hydrogen as + , and empty lattice sites 
as •. 

associated phase transitions give rise to sharp features in experimentally mea- 
sured adsorption isotherms and voltammograms. The static aspects of these 
experimental features are well reproduced by our lattice-gas models with ef- 
fective interaction constants. The interactions in the physical system arise 
from a number of different sources, including mediation through the substrate 
electrons, through phonons, and through the fluid near the surface. Their 
calculation from hrst principles is therefore still not feasible in general. The 
alternative route followed here provides a microscopic picture of the adsorbate 
structure, as well as a procedure for estimating approximate effective interac- 
tion energies from experimentally observed thermodynamic quantities. These 
effective interactions are generally tahen as the minimal set necessary to en- 
sure the stability of the observed ordered phases. Although one cannot expect 
such estimates to be highly accurate, we find that they lead to models with 
considerable predictive power regarding the dependences of observed thermo- 
dynamic quantities on bulh solute concentrations and, to a somewhat lesser 
degree, on temperature. Our studies further demonstrate the importance of 
the geometric fit between substrate and adsorbate. This provides a power- 
ful means of predicting from considerations of local bonding geometry which 
substrate-adsorbate combinations are most likely to give rise to commensurate 
ordered adsorbate phases. Since our methods are not very computationally in- 
tensive, they should be well suited for data analysis by experimental groups. 

Lattice-gas modeling along the lines described in this paper provides basic 
information about the structure and stability of ordered adsorbate phases on 
single-crystal electrode surfaces, that is not readily available by other the- 
oretical or numerical methods at the present time. Such knowledge in turn 
is necessary to develop techniques to manipulate the adsorption of ordered 
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phases in applied fields, such as electrocatalysis. 



Acknowledgement 

We acknowledge useful discussions with L. Blum, K. Franaszczuk, M. Ito, 
A. Tadjeddine, M.F. Toney, and J.X. Wang, correspondence with K. Itaya, 
comments on the manuscript by B.M. Gorman, CCA. Giinther, M.A. Novotny, 
and S.W. Sides, and technical and artistic assistance with Ref. [19] and Fig. 3 
by E. Pepke and D. Poindexter. We thank the authors of Ref. [43] for allow- 
ing us to use their data before publication. This research was supported by 
the Florida State University Supercomputer Computations Research Institute 
under US Department of Energy Contract No. DE-FC05-85ER25000, by the 
Florida State University Center for Materials Science and Technology, and 
by the University of Illinois Frederick Seitz Materials Research Laboratory 
under US Department of Energy Contract No. DE-AC02-76ER01198. Work 
at FSU was also supported by US National Science Foundation Grants No. 
DMR-9013107 and DMR-9315969. 

References 

[1] See, e.g., A. Zangwill, Physics at Surfaces (Cambridge Univ. Press, Cambridge, 
1988), Ch. 11, and references cited therein. 

[2] L. Blum, Adv. Chem. Phys. 78 (1990) 171. 

[3] P.A. Rikvold, J.B. Collins, G.D. Hansen, and J.D. Gunton, Surf. Sci. 203 (1988) 
500. 

[4] J.B. Collins, P. Sacramento, P.A. Rikvold, and J.D. Gunton, Surf. Sci. 221 
(1989) 277. 

[5] P.A. Rikvold and M.R. Deakin, Surf. Sci. 294 (1991) 180. 
[6] P.A. Rikvold, Electrochim. Acta 36 (1991) 1689. 

[7] D. Armand and M. L. Rosinberg, J. Electroanal. Chem. 302 (1991) 191. 

[8] D.A. Huckaby and L. Blum, J. Electroanal. Chem. 315 (1991) 255. 

[9] L. Blum and D.A. Huckaby, J. Chem. Phys. 94 (1991) 887. 
[10] L. Blum and D.A. Huckaby, J. Electroanal. Chem., in press. 
[11] L. Blum, M. Legault, and P. Turq, J. Electroanal. Chem., in press. 
[12] L. Blum and D.A. Huckaby, in this issue. 



15 



[13] A. A. Kornyshev and I. Vilfan, Forschungszentrum Jiilich preprint (1994), and 
references cited therein. 

[14] M. Rubel, C.K. Rhee, A. Wieckowski, and P.A. Rikvold, J. Electroanal. Chem. 
315 (1991) 301. 

[15] P.A. Rikvold and A. Wieckowski, Phys. Scr. T44 (1992) 71. 

[16] P.A. Rikvold, A. Wieckowski, and Q. Wang, Electrochem. Soc. Conf. Proc. Ser. 
93-5 (1993) 198. 

[17] P.A. Rikvold, A. Wieckowski, Q. Wang, C.K. Rhee, and M. Gamboa, in: 
Computer Simulation Studies in Condensed Matter Physics VI, Eds. D.P. 
Landau, K.K. Mon, and H.B. Schiittler (Springer, Berlin, 1993). p. 162. 

[18] M. Gamboa-Aldeco, P. Mrozek, C.K. Rhee, A. Wieckowski, P.A. Rikvold, and 
Q. Wang, Surf. Sci. Lett. 297 (1993) L135. 

[19] L. Hightower, Q. Wang, P.A. Rikvold, M. Gamboa-Aldeco, and A. Wieckowski, 
Urea Adsorption on Platinum (100) Electrodes: A Monte Carlo Simulation, 
FSU-SCRI Video (1993). Electronically available in MPEG format via ftp 
from ftp.scri.fsu.edu/pub/rikvold/uptlOO.mpg or via World Wide Web at 
http://www.scri.fsu.edu/~vinals/materials_science.html. 

[20] M. Gamboa-Aldeco, C.K. Rhee, A. Nahle, Q. Wang, J. Zhang, H.L. Richards, 
P.A. Rikvold, and A. Wieckowski, Electrochem. Soc. Conf. Proc. Ser., in press. 

[21] J. Zhang et al., in preparation. 

[22] M. Gamboa-Aldeco et al., in preparation. 

[23] E. Domany, M. Schick, J.S. Walker, and R.B. Griffiths, Phys. Rev. B 18 (1978) 
2209. 

[24] E. Domany and M. Schick, Phys. Rev. B 20 (1979) 3828. 
[25] M. Schick, Prog. Surf. Sci. 11 (1981) 245. 

[26] See, e.g., K. Binder, in: Monte Carlo Methods in Statistical Physics, Second 
Edition, Ed. K. Binder (Springer, Berlin, 1986). 

[27] See, e.g., T.L. Einstein, Langmuir 7 (1991) 2520. 

[28] See, e.g., T. Aukrust, M.A. Novotny, P.A. Rikvold, and D.P. Landau, Phys. 
Rev. B 41 (1990) 8772. 

[29] See, e.g., P.A. Rikvold and M.A. Novotny, in: Condensed Matter Theories, 
Vol. 8, Ed. L. Blum and F.B. Malik (Plenum, New York, 1993). 

[30] A. Itai, H. Yamada, T. Okamoto, and Y. Iitaka, Acta Cryst. B33 (1977) 1816. 

[31] See, e.g., Table 4 of C. Kittel, Introduction to Solid State Physics, Sixth Ed. 
(Wiley, New York, 1986). 

[32] A. Tadjeddine and A. Peremans, in this issue. 



16 



[33] D.K. Hilton, B.M. Gorman, P.A. Rikvold, and M.A. Novotny, Phys. Rev. B 46 
(1982) 381. 

[34] See, e.g., M.P. Nightingale, in: Finite-Size Scaling and Numerical Simulation 
of Statistical Systems, Ed. V. Privman (World Scientific, Singapore, 1990). 

[35] W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.P. Flannery, Numerical 
Recipes, Second Ed. (Cambridge Univ. Press, Cambridge, 1992), Ch. 12. 

[36] V. Privman, in: Trends in Statistical Physics (Indian Council for Scientific 
Information, New Delhi, 1994). 

[37] P.A. Rikvold and B.M. Gorman, in: Annual Reviews of Computational Physics 
I, Ed. D. Stauffer (World Scientific, Singapore, 1994). 

[38] M. Acharyya and B.K. Chakrabarti, in: Annual Reviews of Computational 
Physics I, Ed. D. Stauffer (World Scientific, Singapore, 1994). 

[39] B.M. Ocko and J. Wang, Phys. Rev. Lett. 65 (1990) 1466. 

[40] O.M. Magnussen, J. Hagebock, J. Hotlos, and R.J. Behm, Faraday Discuss. 94 
(1992) 329. 

[41] G.J. Edens, X. Gao, and M.J. Weaver, J. Electroanal. Chem., in press. 
[42] K. Itaya, private communication. 

[43] Y.-E. Sung, S. Thomas, and A. Wieckowski, in preparation. 

[44] P. Mrozek, M. Han, Y.-E. Sung, and A. Wieckowski, Surf. Sci., in press. 

[45] C. Pascard-Billy, Acta Cryst. 18 (1965) 829. 

[46] M. Ito, private communication. 

[47] R.J. Baxter, Exactly Solved Models in Statistical Mechanics (Academic, 
London, 1982). 



17 



